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ABSTRACT 

This study employs SeaWiFS data over the waters off the southeastern China to evaluate a semi-analytical algorithm ior 
euphotic zone depth (Z e ). This algorithm is based on water’s inherent optical properties (IOPs), which can be 
near-analytically calculated from spectral remote-sensing reflectance, where remote-sensing reflectance can be derived 
from the normalized water-leaving radiance provided by SeaWiFS. In the Taiwan Strait, compared with in situ Z e (±3 
hour within SeaWiFS collection), average error (e) is 15.0 % and root mean square error (RMSE) is 0.074, with Z e in a 
range of 14-34 m from field measurements. In the South China Sea, compared with in situ Z e (±48 hour within SeaWiFS 
collection),e is 5.1 % in summer and 22.6 in winter, while RMSE is 0.032 in summer and 0.129 in winter, with Z c in a 
range of 10-82 m from field measurements. For comparison, we also evaluate the performance of the empirical Z e 
algorithm that is based on chlorophyll concentration. It is found that the IOP-centered approach has higher accuracy 
compared to the chlorophyll-a centered approach (e.g. in the South China Sea in winter, e is 55.3 % and RMSE is 0.219). 
The new algorithm is thus found not only worked well with waters of the Gulf of Mexico, Monterey Bay and the Arabian 
Sea, but also worked well with waters of the China Sea. 
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Tie euphotic zone is a layer where light is sufficient for phytoplankton production by photosynthesis to exceed the loss 
ofinaterial that takes place through respiration \ Although the disphotic zone has attracted increasing attention recently 2 , 
the euphotic zone is nevertheless the most important in the context of ecosystem dynamics 3 ' 4 and air-sea interaction 
though transfer of either heat 5 or materials, in particular green house gases such as C0 2 6,/ . In practice, the euphotic 
zone depth (Z e ) is defined as the depth at which PAR (photosynthetic available radiation) value is 1% of the surface 
^abe 8 . Similar as Secchi disk depth, Z e is a measure of water clarity and is much more robust. Because Z c is a 
cumulative measure of biogeochemical properties of the upper-water column, changes of Z e also depict environmental 
patterns that might be tightly associated with the climate. 

There are various approaches to estimate Z e from ocean color remote sensing. One simple empirical method is 
based on Case-1 water assumptions 9 ’ 10 , where Z e is calculated from remotely derived concentration of chlorophyll-a 
([Chi]) 9 . Recently, Lee et al. 11 developed an analytical model to describe the vertical attenuation of downwelling vector 
imdiance in the visible domain. Based on this model, Z e can be easily calculated when IOPs (the absorption and 
backscattering coefficients at 490 nm, in particular) are known, either from in situ measurements or from remote sensing 
of ocean color 12 . This IOP-centered approach was evaluated with ship-borne measurements made over three different 
regions (the Arabian Sea, the Monterey Bay and the Gulf of Mexico) at different seasons. It was found that the Z e values 
calculated via Rrs-derived IOPs were within -14% of the measured values on average, while the error was -33% when 
Z e was calculated via Rrs-derived [Chi]. However, how this approach performs in other regions keeps unknown. In this 
study, we test the performances of the IOP-centered approach 12 and the [Chi] approach 9 for waters off the southeastern 
China. More importantly, to the best of our knowledge, it is the first time that both approaches are tested with match-up 
data assembled from in situ measurements and from the SeaWiFS ocean color. The calculated Z c are compared with' 
those from profiles of PAR measurements to evaluate the performance of the approaches. 


2. BACKGROUND OF THE REGION UNDER STUDY 

The region we do the tests for the Z c performance includes the Taiwan Strait and the South China Sea (Fig. 1). The 
Taiwan Strait (TWS), with an average water depth of 60 m, is a shallow shelf-channel that connects the South China Sea 
(SCS) with the East China Sea (ECS). Warm, saline, and oligotrophic water enters the TWS from the SCS while cold, 
fresh, and eutrophic water intrudes into the TWS from the ECS along the Chinese coast. Their relative influence, which 
varies seasonally in response to changes in the monsoonal wind, is a major determining factor of the hydrographic 
conditions and biological productivity in the TWS. Coastal upwelling occurs regularly in summer in this region 13 , 
affecting primary production and fisheries. The near-shore waters receive inputs of coastal runoff, including Minjiang, 
Jinjiang, Jiulong Jiang and Hanjiang rivers with high loads of nutrients and other terrestrial substances. 

The South China Sea (SCS) is one of the major marginal seas. The Pearl River discharges into its northeast, through 
which the SCS receives freshwater as well as nutrients and pollutants from one of the most industrialized regions of 
China. Climatic variations in the atmosphere and the upper ocean of the SCS are primarily controlled by the East Asian 
monsoon, which follows closely the variations in the equatorial Pacific 14 . Recent arguments emerged towards a role of 
C0 2 source the SCS played 15,16 . 
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Figure 1. The map of the region for test of Z e algorithm performance 


3. DATA AND METHODS 


3.1 Z e from in situ PAR 

Field measurements of PAR were carried out in August 1998 and 1999 and July-August 2004 for the TWS, and in 
February 2004 and July and December 2006 for the SCS. Figure 2 shows the sampling stations. Instantaneous PAR 
(400-700 nm) in the upper water column at time t (PARt(z)) and depth z was measured by a PAR sensor (Biospherical 
Instruments, Inc.) mounted on a SeaBird CTD rosette. During these measurements the above-surface solar zenith angles 
(0.) were between 8° and 77° (55 measurements in total, of which 16 (27%) had 0 a > 50°). 

Following Lee et al. 12 , Z e was determined by calculating the ratio of PAR at depth z to surface PAR (PAR(0)) as 

below: 


r 

' PAR(O) 


( 1 ) 


Z e is the depth z where r PAR (z) = 0.01,. Because depth was rarely recorded with r PAR (z) exactly 1%, an approximation was 
obtained by exponentially interpolating r PAR (z) between 0.9% and 1.1%. For some measurements in the SCS, the PAR 
sensor did not provide a reading when PAR(z) was 1% of PAR(0). In that case, only the depth z of 10% of PAR(0) was 
directly obtained from PAR profiles. For the other stations where both depths 10% and 1% of PAR(0) (Z e & Z i0 %) were 
obtained from PAR(z), it was found that a good polynomial relationship existed between Z e and Z )0 % (Fig.3). The 
following function was derived: 

Z e = - 0.0194*Zio% 2 + 2.7619* Z w/ „ - 3.6721 (4.6 <Z 10 %<56.9 m) (2) 
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nriis relationship is slightly different from that of Lee et al. 12 , which could be due to that our data here covers a much 
Avider range of Z 10 o /o (4.6 - 57 m instead of 2.1 - 29 m). Therefore, for those stations where no Z e was obtainable from 
Xpar(z), Z c values were estimated by applying equation (2). 



Figure 2. Location of stations for PAR measurements in the Taiwan Strait (a) and the South China Sea (b). In Fig.2b, the red and the 
blue crosses marked the stations sampled in February and July 2004, respectively. 



Figure 3. Relationship between Z c and Z ]0 o /o from in situ PAR profiles in the South China Sea 

3.2 Remotely sensed Z e 

For the dates having field PAR data, remote sensing reflectance (R rs ) were calculated from SeaWiFS daily Level-2 data 
which were obtained from the NASA Distributed Active Archive Center (DAAC). The resolution was 1 km for the TWS 
and 4 km for the SCS. They were used in the OC4v4 band ratio algorithm 17 to calculate [Chi], and then Z e was derived 
from: 
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Z ( = a([Chl]) b , 


(3) 


where the a value was 34.0 and the b value was -0.39 9 . 

In the meanwhile, the absorption and backscattering coefficients at 490 nm, a(490) and £*(490), were derived from ^ 
R rs by applying the updated QAA 12,18 . Kyis(z), the attenuation coefficient of Eyis (downwelling irradiance in the visible 
domain, 350-700 nm), was then calculated as a function of a(490), £*(490) and the sun angle (0 a ) "■ The IOP-centered 
Z e 12 was thus derived based on the following equation: 

K VIS (Z e )-Z e =4.605 (4 ) 


3.3 Comparison of field data with remote sensed data 

Due to frequent cloud cover in the study region, matching data pairs with in situ data collected within ±2-3h of the 
satellite overpass for a rigorous comparison 19 were limited, especially in the SCS. Matching pairs having the time 
differences of ±3 h at the corresponding cruise survey locations were collected for the TWS. However, for the SCS a 
window of ±48 h was chosen in order to obtain enough matching pairs to ensure statistically meaningful results. 

For simplification, hereafter the Z c values derived from field PAR data were annotated as Z e , those from the 
[Chl]-centered approach as Z e c and from the IOP-centered approach as Z/. 

Following Lee et al. 12 , we used the root mean square in log scale (RMSE) and an averaged percentage error (e) as 
measures to describe the similarity/difference between the in situ and satellite data sets. 


8 = 


-1 

v n >=i 


(4 


(Oi 


x 100% 


(5) 


RMSE = Ji|(log((4 ),) - log«Z"), Df 

RMSE and s were calculated for Z e c by substituting Z/ with Z e c in Equation 5 & 6. 


( 6 ) 


4. RESULTS AND DISCUSSION 


Figures 4-6 showed Z e M versus Z e c and Z/, respectively. Table 1 summarized the error estimates, both in our study region 
and in the Monterey Bay, the Gulf of Mexico and the Arabian Sea (abbreviated as MAG) 12 . 

All the three cruises in the TWS were carried out in summer. The averaged percentage error (s) between Z e M and Zj 
was 15.0% (with a maximum error of 53.6%) for a range of ~ 16—34 m from PAR measurements. By contrast, e between 
Z e M and Z e c was 40.1% while the maximum error was 99.1%. The root-mean-square error in log scale (RMSE) for Zj 
was 0.074, significantly smaller than the RMSE (0.159) of Z e c . 

For the summer case of the SCS (July 2004), the performance of Zj and Z e c was close to each other, e was 5.1% for 
Zj (maximum error 17.1%) while it was 7.2% for Z e c (maximum error 21.3%). RMSE was 0.032 for the former and 
0.041 for the latter. Nevertheless, the winter case of the SCS (Feb.2004) was different, e for Z e c (55.3%) was more than 
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double of that for Zj (22.6%). RMSE was 0.219 for Z e c whereas it was 0.129 for Zj. 

It was obvious that the performance of the IOP-centered approach was better than that of the [Chl]-centered 
approach for data in this study. The average error and RMSE were comparable to those published in the Ref. 12. Our 
region, the southern part of the China Sea, covered a wide range of waters, from very shallow near shore water (- 30 m 
depth) to the deep basin (-2000 m depth). The range of surface [Chi] was -0.01-6.29 mg/m 3 with Z c ranging from -16 to 
80 m. This result was particularly encouraging because we used SeaWiFS Rrs data to evaluate algorithms that were 
developed independently. 



approach, (b) from [Chi]-centered approach. 




Figure 5. Comparison between satellite Z c and the Z c derived from in situ PAR profile in the South China Sea in July 2004 (a) from 
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IOP-centered approach, (b) from [Chl]-centered approach. 



Figure 6. Comparison between satellite Z c and the Z c derived from in situ PAR profile in the South China Sea in February 2004 (a) 

from IOP-centered approach, (b) from [Chl]-centered approach. 
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Table 1 Error characters for the comparison between the Z e derived from in situ PAR profiles and from SeaWiFS Rrs 




IOP-Centered Approach 


Chi-Centered Approach 


n 


References 



Average 

error(%) 


Max.error 

(%> 


RMSE 


Average 

error(%) 


Max.error 

(%) 


RMSE 

TWS 

(summer) 


15.0 


53.6 


0.074 


40.1 


99.1 


0.159 


2 

0 


This 

study 

SCS 

(summer) 


5.1 


17.1 


0.032 


7.2 


21.3 


0.041 


1 

1 


This 

study 

SCS 

(winter) 


22.6 


43.9 


0.129 


55.3 


153.4 


0.219 


2 

4 


This 

study 

SCS 

(autumn) 


6.6* 


N/A 


N/A 


10.7 


N/A 


N/A 


5 


21 

MAG 


13.7 


63.5 


0.079 


32.7 


218 


0.162 


6 

4 


12 


* An empirical algorithm based on Kd(490) 21 rather than the IOP-centered approach 12 was applied. 


The performance of the IOP-centered approach in the TWS looked not as good as that in the SCS in summer, and it 
also appeared worse in winter than in summer in the SCS. One potential reason for the excellent performance in the SCS 
in summer was that all of the 11 data points (there were 8 measurements at one station to obtain a time-series observation) 
obtained in clear oligotrophic water, where IOPs derived from Rrs could be more accurate and the waters could be more 
uniform. Even for the station relatively near shore (see Fig.2b for the location), [Chi] was about 0.2 mg/m 3 at 75 m 
(generally the depth of the subsurface [Chi] maximum) 20 . On the contrary, all the TWS data were collected in shallow 
and relatively optically complex waters, where there could be more spatial variation and it is harder to achieve accurate 
atmospheric correction. In addition, for 62% of the in situ measurements in the SCS in winter the sun angles were >50° 
while 29% of them were >70°. The IOP-centered approach, however, was developed with sun angles as high as 60°. The 
results here suggest that fine adjustment to the IOP-centered approach might be necessary. Nevertheless, the overall 
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excellent agreement between SeaWiFS Z c and in situ Z e provide initial but significant indication that euphotic zone depth 
of wile areas can be well derived from satellite data (such as SeaWiFS Rrs) with the IOP-centered approach, although 
moxeevaluations with in situ observations are certainly necessary and helpful. 

Tang et al. 21 recently developed an empirical algorithm for Z c based on IQ(490) of the SCS, which calculated Z c 
fromMODIS ocean color data (K<j(490)). They compared the calculated Z e with the result from [Chl]-centered approach 
and the in situ Z e from SPMR for measurements made in September 2004. Due to measurement limitations, the data 
coveied in that study had only five points with Z c in a range of 30-70 m, which may not be enough to show statistical 
significance. 
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